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Abstract. Motivated by algorithmic problems arising in quantum field theories 
whose dynamical variables are geometric in nature, we provide a quantum algorithm 
that efficiently approximates the colored Jones polynomial. The construction is 
based on the complete solution of Chern-Simons topological quantum field theory 
and its connection to Wess-Zumino-Witten conformal field theory. The colored Jones 
polynomial is expressed as the expectation value of the evolution of the (/-deformed 
spin-network quantum automaton. A quantum circuit is constructed capable of 
simulating the automaton and hence of computing such expectation value. The 
latter is efficiently approximated using a standard sampling procedure in quantum 
computation. 
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1. Introduction 

A new frontier of quantum information is the search for algorithms capable of addressing 
problems in low dimensional geometry and topology. The Jones polynomial pQ 
characterizes the topology of knots and links (collections of circles in 3-space) and is 
associated with the expectation value of a Wilson loop operator in quantum Chern- 
Simons field theory in three dimensions. The algebraic content of this theory is encoded 
into a quantum group structure. The Jones polynomial is the link invariant obtained 
with all the component knots labeled with the fundamental irrep of the quantum 
deformation of SU(2), denoted in the following by SU(2) q . Efficient quantum algorithms 
for approximating the Jones polynomial have been recently proposed in [2|. In [3] we 
introduced the g-deformed spin network automaton model. The spin-network quantum 
simulator model, which essentially encodes the (quantum deformed) £77(2) Racah- 
Wigner tensor algebra, was shown |3 ; to be capable of implementing families of finite- 
states and discrete-time quantum automata which accept the language generated by the 
braid group, and whose transition amplitudes are indeed colored Jones polynomials. The 
latter are an extension of the Jones polynomial with arbitrary irreps of SU(2) q labeling 
the component knots. In this paper we shall explicitly construct a quantum circuit 
which efficiently simulates the dynamics of these automata and hence, if appropriately 
sampled with a set of measurements, approximates the colored Jones polynomial. We 
shall discuss the complexity of the circuit showing that, since the time complexity of the 
spin network automaton is polynomial in the size of the input (depending on the index 
of the braid group and on the number of crossings of the knot diagram), the algorithm 
that efficiently simulates the automata also provides an efficient estimation of the link 
invariant. 

The paper is organized as follows. In section 2 we briefly review the physical setting 
of quantum geometry and the role that in it play the topological invariants. In section 3 
we concisely describe the structure of the spin network quantum automaton. In section 
4 we provide the details of the quantum algorithm that approximates the colored Jones 
polynomials and of the corresponding quantum circuit. In section 5 we provide a few 
concluding remarks and discuss possible future developments and extensions of the 
methods and concepts introduced in the paper. 

2. Quantum geometry and topological invariants 

General relativity - the prototype of physical theories whose dynamical variable, the 
gravitational field, is geometric in nature - still represents a major improvement in 
the 'geometrization programme' stated by Klein and Einstein almost one century ago. 
These ideas laid dormant long after the birth of quantum mechanics and quantum field 
theory. In particular, the quest for a quantum gravity theory dates back to the sixties, 
when Arnowitt, Deser and Misner j3] introduced the so-called (3 + 1) decomposition of 
Einstein field equations, a hamiltonian reformulation of general relativity to be assumed 
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as the basic ingredient for constructing a canonical quantization scheme for gravity. 
We refer the reader to the classical textbooks E] for accounts on quantum general 
relativistic theories up to the seventies. 

Nowadays such approach has been almost abandoned in favor of (hopefully) more 
effective quantization schemes, but a number of substantial contributions developed 
in that golden age keep circulating. A good example is provided by Wheeler's 
'geometrodynamics', which embodies the concept of 'quantum geometry' of the physical 
three-dimensional space, to be thought of as quantum fluctuations of (diffeomorphism 
classes of) 3-metrics within the 'superspace' [7j. As we shall see below, 3-dimensional 
extended objects - more precisely, smooth 3-manifolds endowed with riemannian or 
lorentzian metric tensors - with their rich geometric structure play a prominent role in 
Chern-Simons quantum field theories and associated statistical field theories. Moreover, 
models of quantum gravity in three spacetime dimensions represent by themselves very 
useful toy models in view of generalizations to the physically significant 4-dimensional 
case. 

Euclidean quantum field theory is the quantization procedure of a classical field 
theory based on 'functional integration', over the space of quantum fluctuations of the 
physical fields {0}, of exp[— S ({$})/ h], where S({(p}) is the classical action defined 
in the Wick-rotated counterpart of Minkowskii spacetime [S]. This approach to 
quantization can be related to classical statistical field theory, and consequently it 
inherits the language and methods proper of statistical mechanics (partition functions, 
phase transitions, etc.). This latter feature is particularly fruitful if some kind of 
discretization prescription is applied to the classical theory and suitably extended to 
the path integrals which turn out to be interpretable as statistical sums or partition 
functionals. Indeed, the most successful quantization scheme for general relativity, the 
'sum over histories', was proposed by Hawking and Hartle 9J borrowing techniques 
from the euclidean path integral approach mentioned above. Its discretized version, 
simplicial quantum gravity, relies on Regge's discrete reformulation of classical general 
relativity [H] and has been widely addressed in the last two decades (see e.g. [TTj IT2] 
and references therein). 

The geometrization programme referred to at the beginning of this section was in 
some sense rephrased as a 'gauge principle' by Yang and Mills in |13j . Non-Abelian 
gauge theories interacting with matter fields and their quantized counterparts still play 
a central role in the physics of fundamental interactions, while pure Yang-Mills theories 
(classical and quantum) were recognized to encode a number of interesting geometric 
features (see e.g. the reviews [HI EE]). 

Within the class of quantum Yang-Mills theories we focus our attention on 
'topological' quantum field theories (TQFT), formulated in terms of axioms by Atiyah in 
[TH] (see also [El EE])- Such theories - quantized through the path integral prescription 
starting from a classical Yang-Mills action defined on an orientable riemannian D- 
dimensional space (time) - are characterized by gauge invariant partition functions and 
observables (correlation functions) depending only on the global structure of the space on 
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which the theories live. The latter geometric functionals are computable by standard 
techniques in quantum field theory and provide novel representations of 'topological 
invariants' for .D-manifolds (and/or for particular submanifols embedded in the ambient 
space) which are of prime interest both in mathematics and in theoretical physics. All 
of this enlightens a new kind of connection between geometry and quantum physics: in 
TQFT the physical degrees of freedom of spacetime geometry are global and not local. 
Four dimensional Einstein gravity quantized through the euclidean path integral is not 
a TQFT, however the role of quantum 3-geometry is once more enhanced since gravity 
in three spacetime dimensions can be reformulated as a gauge theory closely related to 
the SU(2) Chern-Simons TQFT [El ED]. 

Without entering into technical details on TQFT in general, let us just recall some 
of the basic ingredients of Chern-Simons quantum field theory. 

The classical SU (2) Chern-Simons action for the 3-sphere S 3 (the simplest 
compact, oriented 3-manifold without boundary) is given by 



where A is the connection 1-form with value in the Lie algebra su(2) of the gauge 
group, k is the coupling constant, d is the exterior differential, A is the wedge product of 
differential forms and the trace is taken over Lie algebra indices. The partition function 
of the quantum theory is obtained from the 'path integral' prescription, by integrating 
the exponential of i times the classical action over the space of gauge-invariant flat 
SU(2) connections (the field variables) according to the formal expression 



where the coupling constant k is constrained to be a positive integer by the gauge- 
invariant quantization procedure. The generating functional (J2J), written for a generic 
compact oriented 3-manifold M 3 with dM 3 = 0, is a global invariant, namely it depends 
only on the topological type of M 3 . This is basically due to the feature that the space 
of solutions of quantum CS theory is finite dimensional [19j . 

The gauge-invariant observables in the quantum CS theory are expectation values 
of Wilson line operators associated with oriented knots (or links) embedded in the 
3-manifold (commonly referred to as Wilson 'loop' operators). Knots and links are 
'colored' with irreducible representations (irreps) of the gauge group SU(2), restricted 
to values ranging over the set {0, 1/2, 1, 3/2, . . . , k/2}. Integer k will be related to the 
deformation parameter q in U q (su(2)), the deformed universal enveloping algebra of 



In particular, the Wilson loop operator associated with a knot K carrying a spin-j 
irreducible representation is defined, for a fixed root of unity q, as (the trace of) the 
holonomy of the connection 1-form A evaluated along the closed loop K G S 3 , namely 




(1) 




(2) 



SU(2), with q = exp(^f ). 




(3) 
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where P denotes path ordering. 

For a link L made of a collection of knots {Ki \ l = 1, s}, each labeled by an 
irrep, the expression of the composite Wilson operator reads 

s 

W^ 2 .., s [L;g] = J] WjJJfjjg]. (4) 
i=i 

In the framework of the path integral quantization procedure, expectation values of 
observables are defined as functional averages weighted with the exponential of the 
classical action. In particular, the functional average of the Wilson operator (J3J) is 

^ N = ^l ff ^H7 tw , (5) 

' ' 1 1 f[DA] expS^'f" 4 ) 

where Scs (A) is the CS action for the 3-sphere given in (0) and the generating functional 
in the denominator is usually normalized to 1. It can be shown that this expectation 
value, which essentially^ coincides with the colored Jones polynomial 122 EH]- 
depends only on the isotopy type of the oriented link L and on the set of irreps {ji, j s }. 
The original Jones polynomial ilj is recovered when a spin-| representation is placed 
on each link component. However the colored link invariants are more effective than 
Jones' in detecting knots, as discussed in j2o"] . 

The colored invariants (0) are the basic objects that will be addressed for 
computational purposes in the rest of this paper. The reader interested in an account 
of their construction through the quantum group approach may refer to sect. 
3), where the issue of (unitary) braid group representations is also considered. In 
the following section we shall use yet another kind of approach j57j, which relies on 
the introduction of the boundary Wess-Zumino-Witten conformal field theory into the 
Chern-Simons setting. Such approach provides a particularly useful presentation of the 
colored Jones polynomials as expectation values of unitary braiding operators in WZW 
theory. 

We leave for the concluding remarks at the end of the paper the discussion of 
possible extensions of the quantum algorithm discussed in next session to the other 
hard problems arising in the theory of closed (hyperbolic) 3-manifolds. 



3. The spin network quantum automaton 

In the first subsection of this section we shall briefly review automata theory and define 
basic concepts of formal language theory. Then we describe the model of quantum 
automaton relevant in the present context: the spin network quantum automaton, which 
provides a natural connection between quantum computation and link invariants [3]. 

| These polynomials are actually invariants of 'framed links', see e.g. (2H1E]- The connection between 
£j 1 ...j a [L;q] and the genuine colored Jones polynomial is Jj 1 ...j s (L, q) = {q _3tu ( i )/ 4 /(q 1 / 2 — q -1 / 2 )} 
£j 1 ...j s [L], once suitable normalizations for the unknots have been chosen. Here w(L) is the writhe 
associated with the planar diagram D(L) of the oriented link L, defined as w(L) = J2 p £ (p)- The 
summation runs over the self crossing points of D(L) and e(p) = ±1 according to simple combinatorial 
rules. The writhe is easily evaluated from the link diagram by simple counting arguments. 
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3.1. Automata theory 

The theory of automata and formal languages addresses in a rigorous way the notions 
of computing machines and computational processes. We review first some of the basic 
concepts. 

If A is an alphabet, made of letters, digits or other symbols, and A* denotes the set 
of all finite sequences of words over A, a language C over A is a subset of A*. The length 
of the word w is denoted by \w\ and Wi is its z'th symbol. The concatenation of two 
words u, v G C is denoted simply by uv. In the fifties Noam Chomsky |2E] introduced a 
four-level hierarchy describing formal languages according to their structure (grammar 
and syntax): regular languages, context-free languages, context-sensitive languages and 
recursively enumerable languages. The processing of each language is inherently related 
to a particular computing model (see e.g. j2H] for an account on formal languages). 
Here we are interested in finite-state automata, the machines able to accept regular 
languages. 

A deterministic finite state automaton consists of a finite set of states S, an input 
alphabet A, a transition function F : S x A — » S, an initial state Sj n and a set of 
accepting states S acc C S. The automaton starts in Si n and reads an input word w from 
left to right. At the z— th step, if the automaton reads the symbol w^, then it updates 
its state to s' = F(s,Wi), where s is the state of the automaton reading u>j. One says 
that the word has been accepted if the final state reached after reading w is in S acc . 

In the case of a non-deterministic finite-state automaton, the transition function is 
defined as a map F : S x A — > V(S), where V(S) is the power set of S. After reading a 
particular symbol, the transition can lead to different states, according to some assigned 
probability distribution. 

Generally speaking, quantum finite-state automata are obtained from their 
classical probabilistic counterparts by moving from the notion of (classical) probability, 
associated with transitions, to quantum probability amplitudes. Computation takes 
place inside the computational Hilbert space through unitary matrices. In the present 
context we shall confine our attention to the so-called measure-once quantum automaton 
[30] . The latter is a 5-tuple M = (Q, E, U, |q ) ? !<!/))> where Q is a finite set of quantum 
states, E is a finite input alphabet with an end-marker symbol # and U (E) : Q — > Q 
is the set of transition functions induced by reading E. The probability amplitude for 
the transition from the state |q) to the state |q') upon reading the symbol o G E is 
therefore (q\U(a)\q') . The state |qo) G Q is the initial configuration of the system, and 
|q/) is an accepting final state. For all states and symbols the function U(a) must be 
represented by unitary operators. The end-marker $ is the last symbol of each input 
word and computation terminates after reading it. At the end of the computation the 
configuration of the automaton is measured, if it is in an accepting state then the input 
is accepted, otherwise it is rejected. The probability amplitude for the automaton of 
accepting the string w is given by fa (w) = (q/| U(w) |q ) > U(w) = -Ylwiei: U(wi): for 
w = 'Y[i w i : ( : ' : denotes ordered product, and for w we used the product symbol to 
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denote concatenations). The explicit form of /m(w) defines the language C accepted 
by that particular automaton. If P denotes the projector over the accepting states, the 
probability for the automaton of accepting the string w is given by Pm(w) = \\P \q w )\\ 2 , 
where \q w ) = U(w)\q ). 

3.2. The q-deformed spin network automaton 

In this subsection we review briefly the structure of the q-deformed spin network 
automaton model, first discussed in [3]. This quantum automaton is an extension of 
the spin network model of computation, introduced in [2] and worked out in [32J, 
constructed on the combinatorics of the Racah-Wigner algebra of the quantum group 
SU(2) q . The q-deformed spin network is a model for a quantum automaton capable of 
processing the braid group language. From now on we shall refer to the model simply 
as spin network, subsuming the use of the deformed algebra. 

The spin-network can be seen as a collection of graphs <5 n (V,E) parametrized by 
an integer n that is a measure of the size of the automaton. For fixed n, to each vertex 
v G V of <& n (V, E) is associated the total Hilbert space Sy^ of the ordered tensor product 
of n irreps of SU{2) q (at q root of unity), together with a particular binary coupling 
scheme b of the n angular momenta = l,...,n) elements of the set J. Different 
vertices correspond to different binary coupling schemes and admit a realization in 
terms of unrooted binary trees whose nodes are labeled with SU(2) q irreps. The edges 
e G E of <3 n (V,E) are associated with unitary evolutions connecting vertices (Hilbert 
spaces) belonging to V. A restriction is imposed on the type of allowed elementary 
unitary evolutions for the states in S)^\ they can be either braid-like or recoupling-like. 
The former are associated with a unitary representation of the braiding between two 
adjacent leaves of the binary tree; the latter are associated with reconfigurations of the 
binary coupling structure of the tree. The graph <3 n (V, E) is constructed in such a way 
that two vertices are connected by an edge if and only if there exist a braid-like or a 
recoupling-like unitary evolution mapping a state of the first vertex to a state of the 
second vertex (see fig. [I]). 




Figure 1. A portion of the spin network graph. Unlabeled trees are associated 
to particular binary coupling schemes for the total Hilbert space. Single 
edges correspond to recoupling-like transformations, double edges correspond 
to braid- like transformations. 
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It was shown in [3| that it is possible to construct a finite-state quantum automaton able 
to process the language generated by the braid group B n . Each graphical realization 
of the quantum automaton can be mapped onto a path in (3 n (V,E). The input-word 
to the automaton is an element b G B n and determines the evolution of the automaton 
according to its image U(b), a unitary representation of B n , which constitutes the 
transition rule. The evolution of the automaton is a sequence of allowed moves on 
<S n (V,E), as depicted in fig.0. 




Figure 2. A portion of the spin network graph with a ribbon denoting a 
path on the graph corresponding to a particular evolution of the spin-network 
automaton. 

The main result in [Hj is that the probability amplitude for the automaton evolution 
associated to the unitary representation of a braid, whose closure is a particular link L, 
is equal to the colored Jones polynomial of L. This result is based on the work of R. 
Kaul j^Zj. The details of the construction of the unitary representations of B n will be 
summarized in the next section. 

The connection between the spin-network computational model, the theory of 
quantum automata and link invariants will allow us to provide a quantum algorithm 
for the efficient approximation of topological invariants of knots. The advantage of 
using the spin network quantum automaton resides in its transition rules, which can 
be straightforwardly expressed in the g-deformed co-algebra recoupling scheme. The 
latter naturally provides the set of unitary operations which are the building blocks 
of the quantum circuit evaluating the invariants. Previous discussion was aimed to 
mapping the problem of evaluating link invariants into the problem of simulating the 
corresponding evolution of the quantum automaton and considering henceforth the two 
problems as equivalent. 

4. A quantum algorithm that approximates the colored Jones polynomial 

In this section we provide a quantum algorithm that efficiently approximates the value of 
the colored Jones polynomial. The interest in this problem stems from the fact that an 
additive approximation of the Jones polynomial is sufficient to simulate any polynomial 
quantum computation The construction of the algorithm involves three different 

contexts: 

(i) a topological context, where the problem is well defined and which allows us to 
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recast the initial instance from the topological language of knot theory to the 
algebraic language of braid group theory; 

(ii) a field theoretic context, where tools from CS topological field theory and WZW 
conformal field theory are used to provide a unitary representation of the braid 
group; 

(iii) a quantum information context, where the basic features of quantum computation 
are used to efficiently solve the original problem formulated in a field theoretic 
language. 

We shall not discuss the topological context itself, where theorems and algorithms are 
available to relate links and braids, and refer the interested reader to [SHI EH] and 
[2E]- The field theoretic context will be discussed in the first subsection. The second 
subsection will deal with the basic structure of the algorithm and its computational 
complexity. In the last subsection we shall complete the proof of efficiency and we shall 
provide notions needed to completely characterize the algorithm. 

4-1- The Kaul construction 

In [21] R. Kaul provides a unitary representation of the braid group and develops a 
method to evaluate observables in SU{2) q CS field theory on the 3-sphere S 3 . His 
construction is based on the relationship between CS theory on a three-manifold 
with boundary and the induced WZW conformal field theory on the boundary. Let 
us consider a three-manifold M 3 with a number n of two dimensional boundaries 
S 1 , S 2 , S n . For each of these boundaries, say S J , there are a number of Wilson lines 
carrying spins intersecting the boundary at some "'puncture'" P { % on the boundary 
(see fig.0). 





Figure 3. Three Wilson lines intersecting the boundaries of a three-sphere. 

We can associate to each S* an Hilbert space Sf . The CS functional integral over 
M 3 is then given as a state in the tensor product of such Hilbert spaces. Following 
the literature, in this section we shall henceforth denote by SU(2)k the quantum 
group SU{2) q with q = exp(|^|); in the following we shall use both expressions 
interchangeably. The conformal blocks of SU(2)k WZW field theory on the boundaries 
S l with punctures determine the properties of Sj l . For each Sj l there are different bases 
related by duality of the correlators of the WZW conformal field theory. 
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Figure 4. Duality transformation between two types of conformal blocks. 



These duality matrices can themselves be expressed in terms of g-deformed SU(2) 
recoupling coefficients. This allows us to relate the unitaries generating the 
computational dynamics of the spin-network automaton to the observables of the CS 
field theory. 

An important aspect in the construction developed in [21] is the close connection 
between links and braids. One obtains this important result by two main theorems. 
The first generalizes to colored oriented braids a theorem, due to Birman [33], relating 
links to plats of braids. The second, which allows us to decompose the duality matrix 
associated to a general q-3nj recoupling transformation into a sequence of elementary 
duality matrices associated to q-6j recoupling transformations, reads: 

Theorem. The correlators for 2m primary fields with spins Ji, J2, J2m in SU(2)k 
Wess-Zumino- Witten conformal field theory on S 2 are related to each other by 



I* 



(p;r)(jl,-,j2m)) = Yl A< S' S r) 
(q;s) 



Jl J2 

33 k 

Hm-\ 32m 



l$(q;S) {jl,-,hm)}, (6) 



where the duality matrix is given as a product of the basic duality coefficients for the 
four-point correlators as 



4 (q;s) 



Jl 32 



ra-1 



m-2 

X 

2=0 



,...,tm- 

1-2 





4r 




) 






Si J2m 





<Jl+l 



U 321+2 
321+3 h+l 



(J) 



Here r = Po,r m - 2 = Pm-i,t = ji,t m -i = hm, «o = ?o, s m -2 = q m -i,hm = 
EilTi" 1 3i an d the spins meeting at trivalent vertices in fig. [3] satisfy the fusion rules of 
the SU(2)k CFT. In fig. [S] we provide a pictorial example of the content of the theorem. 
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q-6j 



q-6j 



q-6j 



q-6j 



|q-6j 



|q-6j 



Figure 5. The sequence of decompositions into elementary recoupling 
transformations of a particular duality matrix. 



The elements in the string ■■■,j m ,jm} will be referred to as j-type numbers, the 

elements in {p , ...,p m _i} as p-type numbers and the elements in {r , r m _ 2 } as r-type 
numbers. 

A general n-strand colored oriented braid is specified by giving n assignments 
ji — Ui, e i), representing the spin and the orientation at each point on the upper and 
lower horizontal lines intersecting the strands. The generators of the groupoid of colored 
oriented braids are 



ji ji+i 



Jl ■■■ Jl+l Jl ■■■ Jn 
jl ••• jl jl+l ••• jn 



(8) 



with / G {!,..., n— 1}, where the "*" implies opposite orientation of the strand with 
respect to the horizontal line (fig. jU]). 



A* A* A* A* A* A* 
J1 J2 jl + 1 jl jn-1 jn 



A A A A A A 
J1 J2 Jl Jl+1 Jn-1 j„ 



Figure 6. A graphic realization of a generator of the colored braid group. 

The generators of colored oriented braids satisfy the usual defining relations of the braid 
group (fig. [7|) 



Ji+l Ji 
Ji Ji+l 



+i 



Ji+2 Ji 
Ji Ji+2 



Ji+2 Ji+l 
ji+l ji+2 
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\ Ji+1 Ji+2 J \ Ji Ji+2 J \ Ji Jt+1 

for % = 1, n — 1 and 

ji ji+i J \ ji ji+i J \ ji ji+i J \ it ji+i 

for \i-l\> 2. 



A* A* A* A* A* A* 

Ji+2 Ji+1 Ji Ji+2 Ji+1 Ji 




AAA AAA 
Ji Ji+1 Ji+2 Ji Ji+1 Ji+2 



Figure 7. Defining relations for the colored braid group generators. 



The platting of a colored oriented braid on an even number of strands is the pairwise 
joining of contiguous strands, both from above and below. Birman's theorem, which 
relates oriented links to plats of ordinary braids is extended in j2Z| to colored 
oriented braids in such a way that a colored oriented link is represented by the plat 



closure of an oriented colored braid b 



n 

3x 



I in 
3m 



r 

rn 

7* 



see fig.jBj. 



A A A A 

I I* I I* 

1 1 2 2 



A A 
I I* 



2m-braid 



platting 



A A 



A 



A 



2m-plat 



A A A A 



J J 



v7\7 \JV 



Figure 8. Platting of 2m colored strands. 

We can finally describe now a method for evaluating the expectation value of an arbitrary 
Wilson link operator. Consider the three-sphere S 3 with two three-balls removed. This 
is a manifold with two boundaries with the topology of the 2-sphere S 2 . Let us place in 
this manifold 2m Wilson lines with spins ji,j2, J2m, such that all the spins generate 
an SU{2) q singlet connecting one boundary to the other. It is easily recognized that 
with these Wilson lines we can realize any element of i?2m (see fig.jH])- 
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Figure 9. An arbitrary colored braid pattern embedded into a three-manifold 
with boundaries, £ , £ 2 , with the topology of S 2 . 



The CS functional integral over the three-manifold can be realized by a state in the 
tensor product of vector spaces f) 1 (g> i} 2 , associated with the two boundaries S 1 and 
X 2 . Conformal blocks can be chosen as basis vectors for these vector spaces. The inner 
products of these basis vectors are normalized according to 

($(P;r)(3T>3'2> -> J2m)l $ (U;V)(jl, 3*2, - J*2m)) = ^p,U^r,V- (9) 

The basis vectors |$(p ; r)(jij2---j2m)) of the conformal blocks $(p ; r)(jij2---j2m) are 
eigenfunctions of the odd indexed braiding generators &2Z+1 defined in (jHJ). The even 
indexed braid generators b 2 i are diagonalized in the basis |$(q ; s)(j'i> •••> J2m))- The 
following eigenvalue equations hold 

62J+i|$(p ; r)(j2Z+i, J2Z+2)) = A P; O2Z+1 , J2Z+2 ) I ^(p ; r) (J2Z+2 , J2Z+1 ) ) , (10) 

M$(q;S)(j2Z, j2Z+l)) = A ?i (j2Z,j2Z+l)|$(q;S)(j2Z+l,j2z))- (11) 

Here |$(p ; r)(jz, jz+i)) = |$(p ; r)0i, ■■■,3i,ji+i, ■■;hm))- The eigenvalues of the braiding 
matrices depend on the relative orientation of the strands and, for right-handed half 
twists (i.e. over-crossings) their value is 



X t (j,t\ = (-) j+i "' g( c J+ c d/2+c min(lii) 



-c t /2 



(12) 



for parallel oriented strands, and 

X t (j,i) =(_)l^|-' g -h- Cl |/2+ Ct /2 ; 



(13) 

if the orientation is anti-parallel. Here Cj is the quadratic Casimir operator equal to 
j (j + 1) for the spin j representation. The eigenvalues ()12|) and f)13|) derive from 
the monodromy properties of the conformal blocks of the corresponding CFT. The 
associated unitary representation of the braid group is provided by the following 
theorem: 

Theorem. A class of representations K : B n — > U(d) from the generators of the 
groupoid of colored oriented braids into the unitary d x d matrices (d = d(n, \b\) ) in the 
basis |$(p ; r)), is given by 

(P';r') 



K 



52Z+1 



121+2 
321+1 



Hl+\ 
321+2 



(P;r) 



A Pi [321+1,321+2 ) <£' , 



(14) 
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and by 



K 



hi 



tei+i J21 
hi 321+1 



(p' ; r') 
(P;r) 



E 

(q;s) 



,4 



(q;s) 
(Pit) 



J2Z-1 J2Z+1 

J2i J2/+2 



Ag ; (j2i, J2Z+1 ) -4(q. s) 



321-1 321 
321+1 321+2 



(15) 



The proof that the defining relations for the braid generators are indeed satisfied can 
be found in This result can be used to prove the next 



Theorem. The expectation value of a Wilson loop operator for an arbitrary link 

it 



1 1* 1 

L presented as a plat closure of a colored oriented braid b | ~ 
generated by a word given in terms of the braid generators, is given by 

m 

V[L;y,q] = n^ + i] 



m 

jm Jr. 



8=1 

x <$ (0;0) (/i,...,O|K 



h 

31 



'in 

7* 

Jm 



l*(0;0)(jl> ->&)>> (16) 



where the multi-index (0; 0) denotes the case in which all the elements in the set of p 
and r -type numbers are equal to 0, while 



x 



qX/2 — q x/2 

qi/2 _ q -i/2 



is the standard notation for the quantum integer. The latter theorem gives us the 
explicit evaluation of the colored polynomial. It can be shown that the Jones polynomial 
corresponds to a spin-| representation living on all the components of the link. 

4-2. The qubit representation 

In this section we show how to efficiently implement on a qubit-register the Kaul unitary 
representation K of the colored braid group. We prove that each unitary matrix of 
K (B 2 m), interpreted as a gate acting on a qubit-register, can be efficiently decomposed 
into a set of universal elementary gates. This is done first encoding into a qubit-register 
the basis vectors used in K and then showing how K (bi) can be efficiently compiled for 
every b { G B 2m . 

Each vector in the basis set {|$(p ; r) (Ji, j*, jm,jm))fi corresponding to the 
conformal block shown in fig.^i], is completely characterized by three sets of quantum 
numbers, p, r and j, fully labeling the irreps of SU(2) q . Recall that the p-, r- and 
j-type numbers belong to the set {0, |,...,|}, where k is the Chern-Simons coupling 
constant. This means that each type of number can be specified using [log 2 {k + 1)] 
qubits, where \r\ denotes the least integer > r. An element of the basis can then be 
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encoded using (4m — 3) x [log 2 (k + 1)] qubits. The register we need to use has to code 
only for the p-type and r-type numbers, implying that only (2m — 3) x |~log 2 (k + 1)] 
qubits are sufficient. On the qubit-register we chose the order shown in fig.|10j. 

11 



po 




Pm-1 


n 




frn-3 



Figure 10. Register of qubits for the Kaul representation. 



|po> 

|p m -i> ■ 
|n> 

|rm-3> 

















K(b2i+i) 

















|po> 

| Pm-1 > 

|n> 

|rm-3> 



Figure 11. The gate realization of the odd-indexed braid generators. 



The odd-indexed braid generators are diagonal matrices in the basis of the 
K-module, therefore there is no problem in implementing their action on the 
quantum register (see fig. [H]). The even-indexed braid generators have a less trivial 
representation (see fig.fT2"j). 




Figure 12. The gate realization of the even-indexed braid generators. 



Resorting to the representation in (fi3j) we need to apply two duality matrices, or 
recoupling transformations, in order to explicitly construct the image of these generators 
under K. Each recoupling transformation can in turn be decomposed into a series of 
elementary quantum 6j transformations using ((7|), see e.g. fig. |13j. 
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|po> 
|pi> 

|P2> 
|P3> 

|n> 



q-6j 



q-6j 



q-6j 



1 lh> . 



|fa> 



q-6j 



q-6j 



q-6j 



q-6j 



|qo> 
|qi> 
|q 2 > 

|q 3 > 

|S1> 




Figure 13. The quantum circuit implementing the decomposition of the eight- 
point conformal block in terms of q-6j gates. The corresponding path on the 
spin network graph is shown. 



It follows that the problem of efficiently compiling the general recoupling 
transformation from the eigenspace of odd-indexed braiding operations to the eigenspace 
of even-indexed braiding operations can be mapped into the easier problem of efficiently 
compiling a single q-6j transformation (fig. ^3]). To this end, note that a q-Gj 
transformation, or the corresponding duality matrix, is a unitary transformation from 
states in the conformal block of fig.|15b] into states of the conformal block of fig.|15b]. 



C-q6j 



|m> 



>- 

2>" 
3>" 
4>" 



-|J1> 
-|j 2 > 
"|j 3 > 



L J3 j 4 J 



m> 



Figure 14. Definition of the controlled q-6j transformation. 




Figure 15. The q-6j transformation. 

Each element of the associated unitary matrix is defined in terms of the g-Racah 
coefficients by the following expressions 



|m)(/|j 



A, 



3i 32 
33 H 
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= (_)&+*+*+i4) ^/[ 2m + l] [21 + 1] ( 31 32 1 ) , (17) 




where all the relevant triplets of (2) ? -irreps satisfy the fusion rules of the WZW CFT. 
Recall that an explicit expression for the g-Racah coefficient 27J is 





= A (juj 2 , 1) A (is, 3a, A (ii, H, m) A (j 2 , j 3 , m) x 
z + 1]! {[x - jx - 32 -l]\[x- j 3 - u -l] ] -[x- ji - j 4 - m\\ 

32 - k ~ m V- bi + h + k + J4 - x] \ \ji +j 3 + l + m-x]\ 
x \j 2 +h + l + m-x]\}~ 1 , (18) 

where [x}\ = [x] [x — 1]! with [0]! = 1, and [•] denotes the q-integer. The sum is restricted 
to all allowed values of x such that the quantum integers entering the factorials are non- 
negative and 



A (a, b, c) 



a + b + c} \ [a — b + c] ! [a + b — c 



[a + b + c+l]\ 

Due to the finiteness of the sum, the coefficients ()17j) can thus be efficiently evaluated 
classically for all the 5'L r (2) 9 -irreps. 

For what concerns the action on the qubit-register, elements f|17|) belong to unitary 
matrices of rank 2^ Iog ^ +1 ^, parametrized by the set j of those quantum numbers which 
remain unchanged along the transformation. The crucial fact to notice here is that the 
dimension of these matrices is independent of the size of our problem, given by the index 
of the braid group and the number of crossings. Since there exist efficient methods to 
approximate unitary matrices of given dimension fJS], there exists a sequence of universal 
gates that efficiently approximates every q-6j transformation. The number of elementary 
q-Gj transformations needed to decompose a general q-3nj recoupling transformation is 
2m — 3, linear in the size of the problem. In conclusion, the Kaul representation K 
associated with an arbitrary colored oriented braid can indeed be efficiently compiled 
on a standard quantum computer. The circuit implementing the decomposition of K is 
shown in fig.|16j. 



4-3. The algorithm 

The general structure of the quantum automaton whose dynamical evolution derived in 
[3] is characterized by probability amplitudes whose values corresponds to observables 
of the CS QFT (colored Jones polynomial), can be finally translated into an efficient 
quantum circuit by resorting to a procedure similar to that adopted by Aharonov, Jones 
and Landau in [2]. 

In section 4.2 we have shown that every unitary matrix belonging to the image of 
representation K is efficiently decomposable into a set of universal gates. We now prove 
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|po> - 

IP1> - 
lp 2 > 

|p m -2> - 
|pm-1>- 

|n> . 



|qi> 



- C-q6j 



|tl> 



- C-q6j . 



|t2> 



- C-q6j - 



|tm-2> 



- C-q6j . 



- C-q6j - 



|Sm-3> 



a 

w - 



. C-q6j . 



- C-q6j 



\*> 



_|qo> 
Jqi> 

|q 2 > 

|q m -2> 
|q m -i> 

[S1> 



|rm-3> 



- C-q6j " 



lqm-2> 



|Sm-3> 



Figure 16. The quantum circuit implementing the general duality 
transformation . 



that the information contained in the expectation value V[L; j; q] given by 

($ ( ;0 )('V"Jm)|K 

can be efficiently accessed by a series of measurements. 

To begin with, we recall that the standard procedure used in quantum computation 
to evaluate the expectation value of a unitary relies on a scheme dubbed Hadamard's 
trick. The latter was applied for the first time in j2] dealing just with the problem of 
evaluating the Jones polynomial. We further recall that the notion of approximation 
used in the present context, formalized in j^S], is that of additive approximation, which 
has the following meaning: given a normalized function f(x), where x denotes an 
instance of the problem in the selected coding, we have an additive approximation 
of its value for each instance x if we can associate to f(x) a random variable Z such 
that 

Pr{\f(x)-Z\<5}>3/4, 

for any 5 > 0. The time needed to achieve the approximation must be polynomial in 
the size of the problem and in 5" 1 . The additive characterization of this approximation 
scheme underlies the fact that the interval [Z — 5, Z + 5], which we want to determine, 
is constructed adding ±5 to Z. It also distinguishes this approximation scheme from 
the standard fully polynomial randomized approximation scheme. The normalization 
adopted for the colored Jones polynomial of a link L is provided by the product, over 
all the link components, of the quantum integer related to the dimension of the SU(2) q 
irreps labeling the knots. The problem we are interested in can now be stated as follow. 

Problem: Approximate colored Jones polynomials (Vl). Given a colored braid 
b G -82m of length t, a coloring c, a positive integer k and a real 5 > 0, we want to sample 
from a random variable Z which is an additive approximation of the absolute value of 



7* 

rn 

7* 

Jin 



l $ (0;0)G"l>->3m)) 



Quantum geometry and quantum algorithms 



19 



the colored Jones polynomial of the plat closure of h, evaluated at q = exp (|^|), such 
that the following condition holds true 

Pr (\V(L;j;q)-Z\< 5) > 3/4. 

Here the coloring c denotes the set of all possibly different irreps of SU(2) q labeling 
the component knots of L. 

In the following we shall provide an efficient quantum algorithm for Vl, which solves 
it in 0{po\j(£ , steps. As in 01117] we need the following two lemmas in order to 
prove the efficiency of the algorithm. 

Lemma 1. Given a quantum circuit U of length 0(poly(n)) ; acting on n qubits, and 
given a pure state |$) which can be prepared in time 0(poly(n)), then it is possible to 
sample in 0{po\j(n)) time from two random variables a and b, valued in Z 2 , in such a 
way that < a + ib >= ($>\U\$). 

Lemma 2. For a sufficiently large N, given a set of random variables {ri\i = 1, iV} 
of average value m and square variance v 



Pr 



N 

/•; — III 

1=1 



> 8 < 2 exp(-N5 2 /(4v)). 



The first lemma, which is essentially a reformulation of the Hadamard's trick, can be 
proved as follows. Introduce a single-qubit ancilla 21 and denote by (2» the Hilbert space 
of the qubits acted on by U. Define the unitary C:2l®0— >2l®0 through the action: 

C(|0)®|$}) = |0>® |$), 

C(|l)®|$)) = |1>® (E/|$». 

Initialize then the ancillary qubit in the state 4g (|0) + |1}) G 21 and prepare the system 
in the initial state |$). The action of C maps the initial state into \^>) G 21 ® 

|*> = -^(|0>|$> + |1>(W)). 

The reduced density matrix p % of the ancilla is thus equal to 

p - Tr & mm - -Tr 6 ^ ^ ^ J--^ ^ x 

= i (I 2 + a x Re($\U\$) + ^Im($|C/|$)) , 

where $ denotes the density matrix |$)($| and a x ,a y are Pauli matrices. 

The mean value of a sequence of measurements of a x will approach Re($|[/| < l ) ), 
whereas the mean value of a sequence of measurements of a y will approach Im($|f/|$). 
The second lemma, which is a modified version of the well known Chernoff bound, 
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ensures us that we can approximate these values polynomially in the number N of 
samplings and in the inverse of the precision S^ 1 . 

Summarizing the qubit model for the Kaul representation can be used to efficiently 
compile a unitary representation of the colored braid group, and a sampling procedure 
can then be used to efficiently estimate the value of the colored Jones polynomial. The 
circuit that realizes all these steps is schematically depicted in fig.fTTj. 



|0>+|1> 

\/2 " 



- M((Tx;CTy) 




Figure 17. The circuit realizing the Hadamard's trick for the Kaul 
representation of the braid b. M{a x ;a y ) denotes quantum measurement of 
either a x or a y . 



In conclusion, the sampling lemma tells us that measurements of a x on the first qubit 
will provide the value for Re (V (L, j, q)), while measurements of a y on the first qubit 
will provide the value for Im (V (L,j, q)). 



5. Conclusions 



The g-deformed spin network model provides the natural setting for a quantum 
automaton capable of processing the braid group language. Coding of information in 
the spin network is done in the frame of the coupling scheme associated with the 'co- 
power' A n (SU(2) q ) (iterated co-product) of the network g-algebra. Such parenthesized 
coding lends itself quite naturally to deal with a number of hard combinatorial 
problems, ranging from finite groups word or isomorphism problems |38[ to the 
evaluation of topological invariants. Work is in progress on problems of the former 
type. Focusing mainly on the latter, for consistency with our introductory physical 
setting, we discuss here briefly the role played by the colored link polynomial introduced 
above in 3-dimensional geometric topology jlOlE] i n view of Thurston's 'geometrization 
programme' [12]. Indeed the possible extension of our quantum algorithm to address 
the (classically computationally hard) problems outlined below would represent a major 
breakthrough both in quantum computation and in the theory of closed (hyperbolic) 
3-manifolds. On the physical side, such an achievement would open the possibility of 
'controlling' the quantum algorithmic complexity of three dimensional quantum gravity 
models. It is worth recalling here the alternative point of view introduced by the 
recent attempt by S. Lloyd of unifying quantum mechanics and gravity; where the very 
geometry of space-time is a construct derived from the underlying quantum computation 
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5.1. Reshetikhin-Turaev quantum invariants of 3-manifolds and their quantum 
complexity 

At its foundation, knot theory is a branch of geometric topology, since it allows us 
to explore 3-dimensional spaces by 'knotting' phenomena, namely embedded knots 
'interact' with the topological structure of the ambient 3-manifold M 3 . The content 
of the latter remark is made more stringent by a theorem which asserts that every 
closed connected orientable 3-manifold can be obtained by Dehn 'surgery' along a 
framed link embedded in the 3-sphere S 3 (we refer to [2H] for definitions and proofs). 
Roughly speaking, a tubular neighborhood of each component of the embedded link L, 
represented by S 1 x D 2 (D 2 being the 2-disk), is removed and replaced by D 2 x S 1 in 
a suitable way, generating the new manifold. Formally 

{S 3 , L) — > M 3 L = S 3 \L. (19) 

It can also be shown that equivalent links, namely links which are ambient isotopic, 
give rise to the same type of 3-manifolds (the manifolds obtained by surgeries in the 
3-sphere along equivalent links are homeomorphic) . 

The idea that the Jones polynomial at a root of unity q can be 'amplified' to achieve 
a 3-manifold quantum invariant dates back to Witten, and was further implemented by 
a number of authors ([22] and references therein). Such invariants, which correspond 
to the partition function (j2j) evaluated for a manifold Ai L , are linear sums of Jones 
polynomials of copies of the link with the components replaced by various parallels 
of the original components. The authors of 03] propose to address the problem of 
designing quantum algorithms for Witten invariants by resorting to Temperley-Lieb 
algebra techniques. 

The quantum algorithm for the colored Jones polynomials discussed in section 4, 
allows us to conjecture that the associated colored 3-manifold quantum invariants at a 
fixed root of unity can be actually evaluated in a quite straightforward way. The explicit 
expression of the (Witten-) Reshetikhin-Turaev quantum invariant for a 3-manifold M. 3 L 
to be used for computational purposes was proposed by Kirby and Melvin [22] and reads 

T{M\;q) = aL^2\ji\S jl ... js [L;q], 
j 

where j stands for the collective assignment of colorings to the link components, the 
summation is over all admissible colorings and [j] = n^i^Ji + !]• £ji-j s \L\ o] is 
given in (JHJ and oll = b nL c UL . Here b and c are numbers depending on the integer k 
(b = \J (2/k) sin | and c = exp[— 2iri(k — 2) /8k]), Ul is the number of link components 
and o~l is the signature of the linking matrix of L. The linking matrix Ml of a framed 
link L is a symmetric matrix whose entry (M£)ij for i ^ j is the linking number between 
components i and j of L. The diagonal elements of Ml are defined to be the integers 
that give the framing of the individual components. The linking matrix , defined here 
in combinatorial terms, is related to the topology of M. 3 L because its determinant (if it 
is non-zero) is the order of the first homology group of the manifolds. 
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5.2. The volume conjecture for hyperbolic 3-manifolds 

In the framework of Thurston's geometrization program jl2j (all three dimensional 
manifolds can be reconstructed starting from eight types of model geometries), 
hyperbolic 3-manifolds play a special role. Recall that a 3-manifold is hyperbolic 
if it is endowed with a complete riemannian metric with constant negative sectional 
curvature. The most interesting case is that of complete hyperbolic manifolds with 
finite volume (the volume being evaluated in the given metric) since the Mostow rigidity 
theorem asserts that any two such manifolds are homeomorphic if and only if they are 
homotopically equivalent (we refer to [45 for details and original references). Typical 
instances of such manifolds are obtained as quotients of the hyperbolic 3-space H 3 by 
discrete subgroups of the full isometry group of H 3 . 

Consider the set of all complete hyperbolic 3-manifolds with a finite volume. Then 
the set of volumes is totally ordered; moreover there exist only finitely many different 
hyperbolic 3-manifolds with the same volume jlH]. Thus the volume of an hyperbolic 
manifold (unlike what happens in the euclidean and elliptic cases) can be considered 
as a topological invariant. Computer geometry is the branch of geometric topology 
devoted to the calculation of these invariants: it still exhibits many open problems 
interesting for a quantum-computational approach but it is most intriguing that such 
a hard 'computational' approach raised discussions among mathematicians about the 
"philosophical" question of the effectiveness of mathematical proofs. 

At first sight the above remarks on hyperbolic volumes does not seem related to 
our central issue of colored quantum invariant of 3-manifolds. However this is not the 
case: a connection can be easily recognized by observing that most manifolds obtained 
by surgery on framed knots (links) in the 3-sphere can be endowed with hyperbolic 
metrics. Let us focus for example on 'hyperbolic knots', namely those knots which 
give rise by surgery to (finite volume) hyperbolic 3-manifolds: the 'volume conjecture' 
proposed by [HI HHI (see also the review [IU] for extended versions) can be cast in this 
case in the form 



where K is a hyperbolic knot and the notation J^(K) stands for the iV-colored 
polynomial of K evaluated at q = exp(2iri/N). 

Notice that all the quantum algorithms dealing with link polynomials are 
established for a fixed choice of the root of unity q appearing in the argument of the 
invariants, while the volume conjecture involves the analysis of the asymptotic behavior 
of 'single-colored' polynomial of the same knot for increasing values of the coloring itself. 
Recently, Aharonov and Arad [49] have addressed an asymptotic analysis (k — > oo) for 
the Jones polynomial, still ^-colored. It would be interesting to explore the possibility of 
borrowing some of their techniques to test the conjecture (|20jl within the computational 
framework designed for colored polynomials. 
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We are in debt with Romesh Kaul for clarifying remarks on his work on colored 
polynomials. 
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